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1.  INTRODUCTION 


This  first  annual  report  describes  work  done  during  the  period  May  6,  1993  through  May 
5,  1994.  The  objective  of  this  program  is  to  develop  a  phenomenological  chemical-kinetics 
database  for  the  nuclear-disturbed  atmosphere.  To  this  end,  the  comparison  of  data  from  the 
EXCEDE  III  experiment  to  a  robust  chemical  kinetics  model  which  contains  the  current  state 
of  knowledge  of  electron-atmospheric  interactions  will  ultimately  provide  the  information 
necessary  for  predicting  nuclear  effects  on  systems.  Current  uncertainties,  in  the 
phenomenological  models  based  on  reaction  rate  coefficients  and  calculated  electron 
concentrations  that  affect  infrared  emissions  are  thought  to  lie  within  a  factor  of  2  or  3  in  the 
best  known  situations  and  much  worse  in  many  other  cases.  Unknowns,  or  disputed  reaction 
rates,  include  fractional  yields  of  species  such  as  N(^P)  and  N2(A),  effects  of  these  species  in  the 
disturbed  atmosphere,  and  creation  mechanisms  for  species  such  as  rotationally  excited  NO. 
Emissions  due  to  very  high  rotation  levels  which  are  not  quenched  in  the  atmosphere  may  have 
a  great  impact  on  the  atmospheric  bandpasses  where  current  surveillance  instruments  are  being 
designed. 

The  major  effort  during  this  period  has  been  the  analysis  of  electron-beam  enhanced 
NO(Av=l)  emission  observed  by  the  three  circular-variable  filter  (CVF)  spectrometers  and 
interferometer  during  EXCEDE  III.  In  addition  to  considering  the  traditional  N(^D)-l-02  reaction 
at  thermal  energies,  the  chemistry  of  translationally  hot  N  atoms  in  NO  formation  has  been 
investigated  in  detail,  with  emphasis  on  the  hyperthermal  N(‘'S)4-02  reaction.  It  has  been  shown 
that  a  quantitative  description  of  rotationally  thermal  and  hot  NO  emission  observed  in  EXCEDE 
III  requires  the  inclusion  of  hyperthermal  N  atom  chemistry  in  our  chemical  kinetics  model. 

Section  2  discusses  the  evidence  for  the  importance  of  the  reactions  of  hyperthermal  N(‘’S) 
and  N(^D)  with  O2  in  producing  NO.  The  chemical  dynamics  of  the  N('‘S)-l-02  reaction  is 
discussed  in  Section  3,  and  preliminary  calculations  of  NO  vibrational  populations  using  a 
hyperthermal  nitrogen  atom  kinetics  model  is  presented  in  Section  4.  A  brief  summary  is  given 
in  Section  5. 
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2.  IMPORTANCE  OF  TRANSLATIONALLY  HOT  NITROGEN  ATOMS 


The  reactions  of  metastable  nitrogen  atoms  with  molecular  oxygen  are  thought  to  be  the 
major  source  of  nitric  oxide  chemiluminescence  in  the  thermosphere.  Earlier  models  of  NO 
formation  assumed  that  the  chemistry  of  N(^D)  is  the  most  significant  contributor  to  NO 
formation,  while  the  relaxation  of  N(^P)  to  NC^D)  and  the  reaction  of  N(^P)  with  O2®  make  a 
relatively  minor  contribution  in  the  thermospheric  dayglow®  and  aurora. 

Analysis  of  the  EXCEDE  III  CVF  NO  data  at  apogee  (115  km)  has  provided  strong 
evidence  for  a  mechanism(s)  of  NO  formation  in  addition  to  the  thermal  reaction 

N("D)  +  O2  ->  NO(v)  +  O  .  (1) 

Reaction  (1)  has  previously  been  considered  as  a  source  of  rotationally  thermal  NO.  The  source 
of  rotationally  hot  NO  is  currently  unknown,  although  it  has  been  suggested  that  the  reaction  of 
N(^P)  with  O2  is  a  possible  source.  We  have  shown  previously  that  this  process 

N(2p)  +  O2  -  NO(v,T,)  +  O  (2) 

can  not  be  a  source  of  rotationally  thermal  or  hot  NO  emission  in  EXCEDE  III.  To  quantify 
this  discussion,  a  comparison  of  the  CVF  data  at  115  km  (apogee)  with  the  EXCEDE  kinetic 
model  using  Reaction  (1)  as  the  only  source  of  NO(Av  =  l)  emission  is  shown  in  Figure  1. 
Although  the  disagreement  between  the  model  and  data  at  0.2  s  (8"  CVF)  and  0.5  s  (20°  CVF) 
could  be  due  to  uncertainties  in  the  atmospheric  profiles,  it  is  clear  that  the  disagreement  at  0 
s  (0°  CVF)  results  from  the  inability  of  Reaction  (1)  to  explain  NO  emission.  Further  evidence 
of  the  inadequacy  of  Reaction  (1)  is  provided  by  the  analysis  of  the  interferometer  data  at  apogee 
which  shows  significant  emission  from  highly  rotationally  excited  NO.  Furthermore,  at  apogee. 
Reaction  (1)  does  not  produce  the  correct  NO  vibrational  distribution  to  explain  the  rotationally 
thermal  populations  obtained  from  the  interferometer  data. 
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Figure  1.  Comparison  of  the  CVF  Data  in  the 
NO(Av=l)  Bandpass  with  the  EXCEDE 
Kinetic  Model  at  1 15  km  Using  the  Thermal 
N(^D)+02  Mechanism. 


Thus,  it  would  appear  that  another  mechanism  for  the  creation  of  NO(v)  is  necessary  to 
explain  the  data  near  apogee.  Recently,  it  has  been  suggested  that  the  reaction  of  translationally 
hot  N(‘*S)  with  O2  is  an  important  contributor  to  NO  formed  in  the  thermosphere.'®'^ 
Furthermore,  Sharma  et  al.'*^  have  performed  a  calculation  that  shows  the  importance  of  the 
N(‘‘S)+02  reaction  in  recent  observations  of  highly  rotationally  excited  NO  vibrational  emissions 
in  the  dayglow.'’  '”’  Their  calculation  was  based  on  two  recent  studies;  a  quasiclassical  trajectory 
study  of  N(^S)  +  O2  by  Gilibert  et  al.'"’  indicating  that  the  N('*S)  reaction  produces  NO  in  high 
ro-vibrational  states  at  a  nearly  gas  kinetic  rate,  and  the  calculated  nonthermal  thermospheric 
N(^S)  translational  energy  distribution  of  Shematovich  et  al."^^ 
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3.  CLASSICAL  DYNAMICS  OF  THE  NCS)+02(X^E-g)-N0(X2n)+0(^P)  REACTION 


3.1  Introduction 

The  only  published  experimental  data  on  the  N(^S)+02  reaction  are  thermal  rate  constant 
measurements/*^^  which  are  only  reliable  up  to  1500  K,  and  the  NO  vibrational  distribution  at 
room  temperature.^*^***  The  dependence  of  the  reaction  rate  and  final  state  distributions  on  the 
initial  translational,  vibrational,  or  rotational  energy  has  only  been  calculated  on  the  ground  state 
potential  energy  surface.  Furthermore,  the  internal  energy  distribution  of  the  NO(X^n)  product 
is  not  completely  characterized  for  hyperthermal  collisional  energies.  The  present  theoretical 
study  determines  the  reaction  rate  constant  and  the  nascent  NO  vibrational-rotational  distributions 
formed  by  the  translationally  hot  N(4S)-l-02  reaction  as  a  function  of  N('*S)  translational  energy 
from  0.5  to  2.75  eV.  The  results  of  this  study  provide  the  reaction  attributes  of  N(‘*S)-l-02 
useful  for  modeling  the  production  of  NO  from  "hot"  N('‘S)^  *’  ***  and  the  NO  dayglow  emission 
spectra.  <*■*”* 

3.2  Previous  Theoretical  Work 

The  reaction  of  N('‘S)-I-02(X  to  form  NO(X^n)-l-0(^P)  is  exothermic  by  1.38  eV  and 

occurs  on  two  electronic  potential  energy  surfaces  of  ^A'  and  ^A'  symmetry  (neglecting 
spin-orbit  coupling).  The  electronic  degeneracies  are  such  that  4/12  of  the  N(‘*S)-l-02  collisions 
occur  on  the  ^A^  surface  compared  to  2/12  on  the  ground  state  ^A'  surface.  The  remaining  6/12 
of  the  collisions  occur  on  the  *A'  surface  which  does  not  correlate  with  the  ground  state  products 
N0(X^n)-l-0(*P),  and  is  not  further  considered.  An  extensive  set  of  ab  initio  quantum 
mechanical  calculations  for  the  ^A'  and  “A'  potential  energy  surfaces  have  been  carried  out  by 
Walch  and  Jaffe**’*  to  characterize  the  saddle  point  geometries  and  to  a  lesser  extent  the 
minimum  energy  path.  The  lowest  ^A'  potential  energy  surface,  which  has  an  estimated 
activation  energy  of  approximately  0.3  eV,**^*  makes  the  dominant  contribution  to  the  reaction 
rate  for  temperatures  below  15{X)  K.  Although  the  ‘*A'  surface,  which  has  an  estimated  barrier 
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of  0.65  does  not  contribute  to  the  reaction  rate  for  temperatures  below  1500  K,  it  must 

be  considered  for  the  translationally  "hot"  N(^S)  reaction  channel. 

There  have  been  two  previous  quasiclassical  trajectory  calculations  reported  for  the 
N(^S)+02  reaction.  The  first  study  was  performed  by  Jaffe  et  using  the  ^A'  and  ‘‘A' 
surfaces  mentioned  above.  However,  their  study  emphasizes  the  thermal  reaction  rate  constant 
at  high  temperatures  and,  although  consistent  with  experiment,  do  not  provide  rate  constants  as 
a  function  of  collision  energy  or  final  state  distributions  of  the  products.  The  second 
quasiclassical  trajectory  study  by  Gilibert  et  al.^"'  only  considered  the  lowest  ^A'  potential  energy 
surface,  and  therefore  is  applicable  to  thermal  collisions  below  1500  K  or  collision  energies  less 
than  approximately  0.75  eV.  More  importantly,  the  fitting  of  an  analytical  functional  form  to 
the  calculated  ab  initio  ^A'  potential  energy  surface  introduced  an  artificial  barrier  of 
approximately  0.55  eV  in  the  reactant  channel,  which  implies  that  the  calculated  threshold  for 
reaction  is  too  high  by  approximately  0.22  eV.  Thermal  rate  constants  calculated  using  the  ^A' 
PES  of  Gilibert  et  al.  would  result  in  an  activation  energy  approximately  5  kcal/mole  higher  than 
experiment.  Therefore,  we  have  carried  out  extensive  quasiclassical  trajectory  calculations^*^ 
using  a  new  analytical  fit  to  the  ^A'  and  ‘*A'  ab  initio  potential  energy  surfaces  of  Walch  and 
Jaffe^*’*  to  predict  the  reaction  attributes  of  N(‘*S)+02. 

3.3  Methodology 

The  most  recent  ab  initio  calculations  of  Walch  and  Jaffe**’*  for  the  ^A'  and  ‘‘A'  potential 
energy  surfaces  (PES)  are  used  as  a  basis  for  the  analytical  representation  of  the  N('‘S)+02 
reaction.  The  study  of  Walch  and  Jaffe  involved  a  total  of  66  ab  initio  points  for  the  ^A'  surface 
and  23  for  the  '‘A'  surface  at  the  complete  active  space  self-consistent  field  (CASSCF)  level 
followed  by  multireference  contracted  configuration  interaction  (CCI)  calculations.  The  "best" 
theoretical  estimate  of  the  barrier  height  for  the  ^A'  surface  was  estimated  to  be  2-3  kcal/mole 
too  high  compared  to  the  experimental  activation  energy.  In  our  model,  the  calculated  ab  initio 
points  in  the  vicinity  of  the  saddle  point  were  adjusted  by  a  constant  factor  such  that  the  barrier 
heights  for  the  ^A'  and  ‘‘A'  are  in  agreement  with  the  experimental  estimate  for  the  ^A'  surface**^* 
and  theoretical  estimate  for  the  '’A'  surface.**’*  The  adjusted  ab  initio  points  were  then  fit  to  an 
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analytical  representation  of  the  potential  energy  surface  suggested  by  Sorbie  and  MurrelF^^  using 
a  nonlinear  least  squares  method.  Care  was  taken  to  ensure  that  the  fitting  process  did  not 
introduce  any  spurious  artifacts  in  the  resulting  analytical  functions.  The  resultant  barrier 
heights  of  the  analytical  functions  are  0.3  eV  and  0.65  eV  for  the  ^A'  and  '‘A'  potential  energy 
surfaces,  respectively.  The  location  and  the  angular  dependence  of  the  barriers  are  also 
accurately  reproduced,  with  a  standard  deviation  of  the  nonlinear  fits  for  both  surfaces  of  0.047 
eV. 

The  quasiclassical  trajectory  (QCT)  method  is  used  to  compute  the  reaction  rate  constants 
and  the  final  vibrational  and  rotational  state  distributions.^*^  The  validity  of  classical  mechanics 
for  calculating  reaction  rate  constants  and  product  distributions  has  been  discussed  elsewhere.^** 
Classical  calculations  are  expected  to  be  reliable  at  energies  greater  than  the  barrier  heights, 
where  quantum  effects  such  as  tunneling  may  be  important,  for  state-to-state  rate  constants  which 
are  strongly  classically  allowed.  Calculations  are  carried  out  separately  for  the  ^A'  and  '‘A' 
potential  energy  surfaces  and  then  averaged  together  with  the  degeneracy  factors  of  ‘/zand  Va, 
respectively.  Standard  Monte  Carlo  techniques  are  used  to  compute  the  thermal  reaction  rate 
constant  as  a  function  of  temperature,  where  the  translational  energy  and  the  initial  (v,  j)  states 
are  chosen  from  Boltzmann  distributions.  Reaction  rate  constants  and  final  (v',  j')  distributions 
are  also  calculated  as  a  function  of  initial  translational  energy  in  the  range  of  0.5-2.75  eV  with 
the  initial  (v,  j)  states  selected  from  a  300  K  Boltzmann  distribution.  The  classical  equations  of 
motion  are  integrated  using  a  variable  step  size  predictor-corrector  method.®^*  The  final 
vibrational  and  rotational  quantum  numbers  (v',  j')  are  obtained  from  the  correspondence  rules 
v'^i  =  (v'  -f-  ‘A)h  and  j  'el  =  (j '  +  where  v'e,  and  j  'ei  are  the  classical  vibrational  and  rotational 

action  variables.  The  final  (v',  J')  distributions  are  then  obtained  using  the  standard  histogram 
method.  A  total  of  580,000  (330,000  for  the  ^A'  PES  and  250,000  for  the  ''A'  PES)  trajectories 
are  used  in  the  current  study. 

3.4  Results  and  Discussion 

Quasiclassical  trajectory  thermal  reaction  rate  constants  are  computed  for  the  temperature 
range  of  500  K  to  5000  K.  A  comparison  of  the  calculated  thermal  rate  constant,  k(T),  with  a 
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Figure  2,  Comparison  of  the  Experimental  (EXP)  and 
Quasiclassical  Trajectory  (QCT)  Thermal 
Rate  Constants  as  a  Function  of 
Temperature. 

recommended  fit  to  the  available  kinetic  data<‘’>  is  shown  in  Figure  2.  The  QCT  thermal  rate 
constant  agrees  well  with  the  experimental  measurements  over  the  entire  temperature  range, 
although  the  calculations  consistently  underpredict  the  rate  constants  below  1250  K.  This 
systematic  disagreement  is  probably  due  to  two  factors:  either  the  barrier  height  of  0.3  eV  is  still 
too  high,  or,  owing  to  the  rather  large  barrier,  there  is  a  significant  contribution  from  quantum 
mechanical  tunneling,  which  cannot  be  accounted  for  with  the  QCT  method.  Although  of  critical 
importance  for  determining  low  temperature  thermal  rate  constants,  neither  one  of  these  concerns 
should  be  as  significant  for  calculation  of  the  high  energy  (>0.5  eV)  rate  constants  and  final 
state  distributions  reported  here.  The  comparison  in  Figure  2  also  points  out  the  difficulty  in 
extrapolating  low  temperature  rate  constants  to  higher  temperature,  where  a  significant 
contribution  to  the  rate  constant  may  come  from  a  mechanism  or  potential  energy  surface  not 
important  at  lower  temperatures.  In  particular,  the  present  calculations  indicate  that  if  the  lower 
temperature  rate  constant  (only  ^A')  is  extrapolated  to  estimate  the  reaction  rate  constant  at  5000 
K,  the  total  rate  constant  would  be  underestimated  by  a  factor  of  2.  It  should  also  be  mentioned 
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ErfeV) 

Figure  3.  Contribution  of  the  QCT  Rate  Constants  for 
the  ^A'  and  '‘A'  Potential  Energy  Surfaces  to 
the  Total  Rate  Constant  as  a  Function  of  the 
Initial  Relative  Translational  Energy. 

that  the  present  results  are  in  good  agreement  with  the  rate  constants  calculated  by  Jaffe  et  al.^°^ 
using  a  completely  different  analytical  form  for  the  potential  energy  surfaces. 

A  more  useful  quantity  for  the  thermospheric  NO  production  models*^  *’  ’*^  is  the  reaction 
rate  constant,  k(Ex),  as  a  function  of  initial  translational  energy,  Ej.  Figure  3  shows  the  total 
reaction  rate  constant  calculated  as  a  function  of  Ej.  The  initial  O2  vibrational  and  rotational 
states  were  averaged  over  a  Boltzmann  distribution  at  300  K.  Also  shown  in  Figure  3  is  the 
individual  contribution  of  the  ^A'  and  '‘A'  potential  energy  surfaces  to  the  rate  constant.  As 
expected,  the  reaction  rate  is  dominated  by  the  ^A'  PES  for  energies  below  1  eV.  As  the  initial 
translational  energy  is  increased,  the  relative  importance  of  the  ^A'  reaction  increases  and 
becomes  greater  than  the  ^A'  reaction  above  1.25  eV  due  to  the  larger  degeneracy  factor. 

In  order  to  assess  the  possibility  of  the  N(‘*S)+02  reaction  contributing  to  the  dayglow  NO 
emission  from  CIRRIS  1A<®  ’>  and  the  NO  emission  from  the  EXCEDE  artificial  auroral 
experiment,^'*^  it  is  important  to  establish  the  NO  vibrational-rotational  distributions.  The 
calculated  NO  vibrational  distribution  at  initial  translational  energies  of  0.5  eV  and  1.5  eV  are 
shown  in  Figure  4.  At  the  lowest  translational  energy  (0.5  eV),  the  vibrational  distribution  is 
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V’ 

Figure  4.  NO  Vibrational  Distribution  as  a  Function  of 
the  Final  Vibrational  Quantum  Number  for 
Initial  Relative  Translational  Energies  of  0.5 
and  1.4  eV. 


relatively  flat  from  v'  =2  to  v'  =6  and  then  rapidly  falls  off  until  the  maximum  vibrational  state 
(v'=9)  is  reached.  As  Ep  is  increased  to  1.5  eV,  the  maximum  vibrational  state  populated  is 
v'  =  14.  The  vibrational  distribution  peaks  around  v'=4  and  then  exponentially  decreases  until 
—  v'  =  12.  Although  the  maximum  vibrational  state  populated  increases  with  energy  and  the 
shape  of  the  distribution  changes,  the  average  vibrational  quantum  number  is  nearly  constant  at 
a  value  of  3.9.  This  behavior  results  because  the  average  of  the  distributions  from  the  two 
potential  energy  surfaces  is  almost  independent  of  energy,  even  though  the  individual 
distributions  show  a  modest  (and  different)  energy  dependence. 

The  present  results  for  the  ^A'  PES  are  in  qualitative  agreement  with  the  calculations  of 
Gilibert  et  al.^"^  The  largest  difference  between  the  calculations  is  our  prediction  that  the 
reaction  threshold  occurs  at  ~0.3  eV,  about  0.2  eV  lower  than  the  previous  calculation.  There 
are  also  small  differences  in  the  NO  final  vibrational  distribution,  as  the  present  calculations 
indicate  that  the  NO  is  formed  with  more  energy  in  vibration.  This  result  is  due  to  our  ^A'  PES 


9 


being  more  "attractive"  than  the  surface  of  Gilibert  et  al.,  which  has  been  shown  to  be  an 
important  factor  in  the  vibrational  excitation  of  the  reaction  products. 

There  have  been  three  previous  measurements  of  the  NO  vibrational  distribution  at  room 
temperature.  The  measurements  of  Rahbee  and  Gibson^’'*^  fall  off  exponentially  with 
increasing  NO  vibrational  quantum  number  in  contrast  with  the  Et=0.5  eV  calculation.  Herm 
et  al.^‘*^  have  pointed  out  the  importance  of  considering  the  effects  of  vibrational  quenching.  The 
experiment  of  Herm  et  al.^’*^  shows  a  NO  vibrational  distribution  which  is  qualitatively  similar 
to  the  QCT  distribution  at  Et=0.5  eV,  although  the  experiment  indicates  a  substantial  population 
in  the  v'  =0  state.  However,  again  due  to  vibrational  relaxation,  it  is  claimed  that  the  measured 
distribution  only  represents  a  lower  bound  to  the  true  nascent  distribution.  An  oscillating  NO 
vibrational  distribution  was  obtained  in  the  most  recent  exijeriment,^'*^  which  disagrees  with 
previous  measurements.  Nevertheless,  the  QCT  calculations  at  0.5  eV  indicate  that  almost  50% 
of  the  available  energy  goes  into  vibration,  somewhat  greater  than  the  experimental  observation 
of  24-34%.  Thus,  the  QCT  calculations  provide  support  for  the  assumption  of  Herm  et  al.^‘*^ 
that  previous  experiments  actually  measured  a  partially  relaxed  vibrational  distribution. 


V’ 

Figure  5.  The  Average  Rotational  Quantum  Number, 
<J'>,  as  a  Function  of  the  Final 
Vibrational  Quantum  Number  v'  at  Initial 
Translational  Energies  of  0.5  and  1.5  eV. 
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Analysis  of  data  from  the  CIRRIS  and  EXCEDE  observations  have  shown 

that  NO  is  formed  in  the  atmosphere  with  extensive  rotational  excitation.  If  the  "hot"  N('‘S)+02 
reaction  is  important  in  NO  formation  as  indicated  by  Sharma  et  al./®^  the  QCT  calculations 
should  show  high  rotational  excitation.  The  calculated  final  average  rotational  quantum  number 
as  a  function  of  final  vibrational  state  is  shown  in  Figure  5  for  translational  energies  of  0.5  eV 
and  1.5  eV.  The  large  amount  of  NO  rotational  excitation  indicated  by  the  calculations  is 
completely  consistent  with  the  field  experiments.  The  NO  rotational  distribution,  summed  over 
all  final  vibrational  states,  can  be  approximated  by  a  Boltzmann  distribution  with  rotational 
temperatures  of  4000  K  and  10000  K  at  translational  energies  of  0.5  eV  and  1.5  eV, 
respectively. 

Finally,  in  the  analysis  of  CIRRIS  lA  and  EXCEDE  III  NO  interferometer  data, 
vibrational  populations  are  determined  assuming  a  single  rotational  temperature  for  the  complete 
vibrational  manifold. Even  though  Figure  5  shows  that  this  is  a  pioor  assumption,  the 
question  remains  concerning  the  sensitivity  of  the  NO  emission  spectrum  to  the  rotational 
distribution.  A  NO  emission  spectrum  at  a  resolution  of  2  cm  '  using  the  vibrational-rotational 
distribution  based  on  the  present  QCT  calculations  is  shown  in  Figure  6.  The  NO 
vibrational-rotational  populations  were  averaged  over  a  uniform  translational  energy  distribution 
from  0.3  eV  to  2  eV.  Ten  rotational  band  heads  from  the  sequence  v'-»  v'-l  with  v'  =  I-^IO  are 
clearly  evident,  in  agreement  with  the  analysis  of  CIRRIS  lA  dayglow  spectra  by  Smith  and 
Ahmadjian.^'”*  A  second  NO  emission  spectrum  was  calculated  using  a  single  effective  rotational 
temperature  of  7000  K  for  each  vibrational  state.  The  agreement  between  this  spectrum  and  the 
one  shown  in  Figure  6  is  excellent,  with  the  largest  difference  of  ~20%  for  the  height  of  the 
band  heads.  Thus,  even  though  each  vibrational  state  possesses  very  different  rotational 
distribution,  the  assumption  of  a  single  rotational  temperature  for  the  complete  vibrational 
manifold  in  the  analysis  of  NO  high  resolution  emission  data  is  well  justified. 
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FREQUENCY  (cm'^) 

Figure  6.  Calculated  NO  Emission  Spectrum  at  a 
Resolution  of  2  cm  '  using  the  Quasiclassical 
Reaction  Rate  for  Each  Vibrational  State 
with  State-Dependent  Rotational 
Temperature. 


4.  PRELIMINARY  ANALYSIS  OF  NO  INTERFEROMETER  DATA 

The  present  trajectory  studies  have  provided  evidence  that  the  process  involving 
hyperthermal  N’^S)  atoms 

N*(''S)  +  O2  NO(v,T,)  -I-  O  (3) 

provides  significant  source  of  NO  which  is  both  highly  rotationally  and  vibrationally  excited. 
In  order  to  access  the  importance  of  Reaction  (3)  in  EXCEDE  III,  we  have  created  a 
time-dependent  kinetics  model  incorporating  a  mechanism  for  producing  NO  from  "hot"  N'C'S). 
The  most  important  inputs  for  this  model  is  the  reaction  cross  section  for  the 
N*(‘‘S)-l-02^N0(v,j)+0  process,  and  the  quenching  of  N*(‘‘S)  by  major  atmospheric  species. 
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The  translational  energy  of  N’(^S)  is  divided  into  6  energy  bins  from  0  to  2  eV  [(>-*0.0625  eV, 
0.0625-*0.125  eV,  0.125-*0.25  eV,  0.25-*0.5  eV,  0.5-*l  eV,  and  l-*2  eV].  The  reaction  rate 
constants  for  N*(‘‘S)4-02  from  the  quasiclassical  trajectory  calculations  are  then  averaged  over 
translational  energy  bins  from  l-*2  eV,  0.5-*l  eV,  and  0.25-*O,5  eV.  The  various  sources,  and 
associated  translational  energy,  of  N*(^S)  from  irradiation  of  the  atmosphere  are 


ep  -H  N2  -  Cp  N(^S)  -H  NeD,^P) 

[0-*2  eV] 

(4) 

-1-  e.  -*  aN("S)  +  (2-a)N(^D) 

[1.7  eV] 

(5) 

O'"  -1-  N2  -*  NO^  +  N(''S) 

[0.75  eV] 

(6) 

N-^  +  O2  -  0^  -H  N("S) 

[1.7  eV] 

(7) 

NO-"  -1-  e,  -*  0  -1-  N(‘‘S) 

[1.5  eV] 

(8) 

0  -1-  N(2D)  -*  0  +  N(^S) 

[1.27  eV] 

(9) 

N+  -h  0  -*  0-^  -f-  N(^S) 

[0.5  eV] 

(10) 

0  -1-  N(2P)  -*0-1-  N('‘S) 

[1.9  eV] 

(11) 

The  translational  energies  from  Reactions  (4)  -  (11)  are  transformed  to  the  N(‘’S)-l-02  center  of 
mass  and  placed  in  the  appropriate  energy  bin.  There  are  estimates  of  the  cross  section  for  the 
quenching  of  N*(^S),  although  the  uncertainties  are  unknown.  We  have  chosen  a  reasonable 
value  of  5xl0  ‘**  cm^  for 

N*("S)  -t-  M  -*  N*(*S)  -I-  M  .  (12) 

which  is  to  be  compared  with  the  gas  kinetic  estimate  of  3xl0  ‘*  cm^.  The  subsequent  cascading 
of  N’('‘S)  translational  energy  from  bin  i  to  i-1  on  successive  collisions  is  simply  scaled  by  the 
relative  velocity  of  each  translational  energy  bin.  It  has  been  assumed  that  each  collision  of 
N*('*S)  with  M  results  in  a  uniform  distribution  of  the  translational  energy. The  rotationally 
hot  vibrational  populations  calculated  from  Reaction  (3)  are  compared  to  the  NO  populations 
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obtained  from  the  interferometer  data^^  at  an  altitude  of  103  km  (upleg)  in  Figure  7.  The 
excellent  agreement  between  the  chemical  kinetics  model  and  data  provide  the  first  quantitative 
evidence  that  Reaction  (3)  is  the  source  of  rotationally  hot  NO  in  disturbed  atmospheres.  Also 
shown  in  Figure  7  is  a  comparison  of  the  vibrational  populations  obtained  from  Reaction  (1) 
using  the  nascent  vibrational  distribution  of  Kennealy  et  al.^*^  with  the  rotationally  thermal 
populations  from  the  interferometer.  Again  the  agreement  of  the  kinetics  model  and  the  data  is 
excellent,  indicating  that  the  N(^D)  is  thermalized  at  103  km  (upleg). 

A  similar  comparison  of  the  rotationally  thermal  and  hot  NO  vibrational  populations  from 
the  model  and  interferometer  data  at  115  km  (apogee)  is  shown  in  Figure  8.  For  the  thermal 
rotational  component,  the  vibrational  populations  from  the  kinetic  model  are  approximately  a 
factor  of  3  below  the  data.  This  comparison  implies  that  the  N(^D)  is  probably  not  thermalized 
at  115  km,  and  the  rate  constant  and  vibrational  distribution  used  for  Reaction  (1)  is  not  correct. 

In  an  attempt  to  account  for  a  hyperthermal  N(^D)+02  contribution  to  NO  formation,  we 
have  created  a  kinetic  model  for  N*(^D)  similar  to  that  used  for  N  (^S).  The  N(^D)+02  rate 
constant  was  assumed  to  have  a  T’^  temperature  dependence,  which  implies  an  energy 
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Figure  7.  Comparison  of  the  NO  Vibrational  Populations 
Derived  from  the  Interferometer  Data  at  an  Altitude 
of  103  km  during  Upleg  with  the  EXCEDE  Kinetic 
Model  Using  the  Nonthermal  N*(^D)+02  and 
N*(‘*S)+02  Mechanisms. 
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independent  cross  section  of  5.7x10  *’  cm’  (i.e.,  appropriate  to  a  reaction  with  no  activation 
energy).  The  vibrational  populations  obtained  from  the  hyperthermal  N(’D)  kinetic  model  are 
compared  with  the  data  in  Figure  9.  It  can  be  seen  that  the  magnitude  of  the  populations  from 
the  model  is  in  better  agreement  with  the  data,  but  the  vibrational  distribution  is  clearly 
incorrect.  However,  recall  that  the  nascent  vibrational  distribution  is  based  on  the  room 
temperature  measurements  of  Kennealy  et  al.*’**,  which  appears  to  be  inappropriate  at  collision 
energies  in  the  1  to  2  eV  range. 

Figure  8  also  contains  a  comparison  of  the  rotationally  hot  vibrational  populations  from 
Reaction  (3)  with  the  interferometer  data.  Although  the  kinetic  model  is  in  reasonable  agreement 
with  the  data,  the  model  underpredicts  the  vibrational  populations  for  low  vibrational  quantum 
numbers.  At  this  altitude  (115  km),  N('‘S)  translational  energies  in  the  range  of  1  to  2  eV  is 
most  impxjrtant  for  producing  NO.  As  is  evident  from  Figure  3,  the  ^A'  PES  makes  a  significant 
contribution  to  the  NO  vibrational  distribution.  Unfortunately,  the  regions  of  the  PES  which  are 
most  important  for  determining  the  vibrational  distribution  on  the  ‘‘A'  PES  are  not  available  from 
the  ab  initio  calculations.  Evidently  the  ‘‘A'  PES  is  more  repulsive  (i.e.,  colder  NO  vibrational 
distribution)  than  the  PES  used  in  Section  3.  Testing  of  this  hypothesis  will  have  to  await  more 
extensive  ab  initio  calculations. 
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Figure  8.  Comparison  of  the  NO  Vibrational  Populations 
Derived  from  the  Interferometer  Data  at  an  Altitude 
of  1 15  km  during  Upleg  with  the  EXCEDE  Kinetic 
Model  Using  the  Thermal  N(^D)+02  and 
Nonthermal  N*(‘*S)+02  Mechanisms. 
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Figure  9.  Comparison  of  the  NO  Vibrational  Populations 
Derived  from  the  Interferometer  Data  at  an  Altitude 
of  1 15  km  during  Upleg  with  the  EXCEDE  Kinetic 
Model  Using  the  Nonthermal  N*(^D)+02  and 
N*(4s)+02  Mechanisms. 
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5.  SUMMARY 


This  report  has  provided  the  first  quantitative  evidence  of  the  importance  of  hyperthermal 
N^S)  and  N(^D)  atoms  in  the  formation  of  vibrationally  and  rotationally  excited  NO.  The  rate 
constants  and  vibrational  distributions  for  rotationally  hot  NO  are  obtained  from  extensive 
quasiclassical  trajectory  calculations  for  the  NCSj+OzCX^Egj-NOCX^nj+O^P)  reaction  using 
realistic  ab  initio  potential  energy  surfaces.  Good  agreement  has  been  obtained  with 
experimental  thermal  rate  constants,  indicating  that  the  calculations  are  a  good  representation  of 
this  reaction.  Reaction  rate  constants  and  internal  NO  vibrational-rotational  distributions  have 
also  been  calculated  as  a  function  of  initial  translational  energy.  The  NO  vibrational-rotational 
distributions  at  low  translational  energy  are  consistent  with  available  experimental  measurements. 
Recent  field  measurement  programs  observed  NO  distributions  with  extensive  vibrational 
excitation  characterized  by  rotational  temperatures  on  the  order  of  4000  K  to  7000  which 
the  present  QCT  calculations  show  are  consistent  with  NO  formation  from  the  N(^S)-t-02 
reaction. 

Quantitative  agreement  between  the  chemical  kinetics  model  developed  for  EXCEDE  and 
the  vibrational  populations  derived  from  the  interferometer  data  is  obtained  under  conditions  of 
thermaliration  of  nitrogen  atoms  (i.e.,  at  103  km  under  max  dose  conditions).  Analysis  of  the 
vibrational  populations  from  the  interferometer  under  other  conditions  indicate  that  hyperthermal 
N(^D)  atoms  also  play  an  important  role  in  NO  formation.  Information  concerning  the 
N(^D)-f-02  reaction  at  hyperthermal  translational  energies,  analogous  to  that  obtained  in  the 
present  work  for  the  N('*S)-l-02  reaction,  is  not  currently  available. 

The  results  of  this  study  indicate  that  a  quantitative  understanding  of  NO  formation  in  the 
thermosphere  will  require  a  detailed  investigation  into  the  dynamics  of  N(^D)-l-02  reaction. 
Limited  ab  initio  studies  of  the  N(^D)-+-02  system  have  been  earned  out  to  characterize  the 
reaction  pathways. Although  a  total  of  five  potential  energy  surfaces  will  have  to  be 
considered  to  treat  the  dynamics  of  N(^D)+02,  the  calculations  are  well  within  the  scope  of  the 
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quasiclassical  trajectory  method  discussed  in  this  report.  Given  a  more  extensive  and  well 
chosen  set  of  ab  initio  points  for  the  PES  than  is  currently  available,  dynamical  calculations  for 
this  system  will  be  pursued  in  this  laboratory  in  the  future. 
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